Vibrational properties of alpha- and sigma-phase Fe-Cr alloy 
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Experimental investigation as well as theoretical calculations, of the Fe-partial phonon density-of- 
states (DOS) for nominally Fe52.5Cr47.5 alloy having (a) a- and (b) a-phase structure were carried 
out. The former at sector 3-ID of the Advanced Photon Source, using the method of nuclear resonant 
inelastic X-ray scattering, and the latter with the direct method [K. Parlinski et al., Phys. Rev. 
Lett. 78, 4063 (1997)]. The characteristic features of phonon DOS, which differentiate one phase 
from the other, were revealed and successfully reproduced by the theory. Various data pertinent 
to the dynamics such as Lamb-Mossbauer factor, /, kinetic energy per atom, E^, and the mean 
force constant, D, were directly derived from the experiment and the theoretical calculations, while 
vibrational specific heat at constant volume, Cy, and vibrational entropy, S were calculated using 
the Fe-partial DOS. Using the values of / and Cv, we determined values for Debye temperatures, 
Qd- An excellent agreement for some quantities derived from experiment and first-principles theory, 
like Cv and quite good one for others like D and S was obtained. 



For many years Fe-Cr alloy system has been of excep- 
tional scientific and technological interest. The former 
stems on one hand from its interesting physical proper- 
ties such as magnetic ones, and on the other hand from 
the fact that it forms a solid solution virithin the whole 
concentration range preserving, at least metastably, the 
same crystallographic structure (bcc). This, in turn, 
gives a unique chance for investigating the influence of 
the composition on various physical properties within 
the same structure as well as adequatly testing differ- 
ent theoretical models and theories. The latter follows 
from the fact that Fe-Cr alloys constitute the basic in- 
gredient of stainless steels (SS) that for a century have 
been one of the most important structural materials [H, 
and, consequently, some properties of SS are inherited 
from the parent alloy. To the latter belongs a tr-phase 
that can precipitate if a quasi-equiatomic Fe-Cr alloy 
undergoes an isothermal annealing in the temperature 
range ~ 800K<T<- 1100 K. The cr-phase has a tetrag- 
onal structure (type D^^ P42/mnm) with 30 atoms dis- 
tributed over five different sites (Table 1). Its physical 
properties are, in general, quite different than those of the 
a-phase of similar composition. Some properties, like the 
magnetic ones, are even dramatically different [2]: other 
properties, like the Debye temperature, seem to be very 
similar The latter is rather unexpected as the hard- 
ness of the CT-phase is by a factor of ~ 3 larger than 



that of the a-phase. To clarify this situation a more de- 
tailed knowledge of vibrational properties of the a- and 
a-phases is essential. In addition, the cr-phase belongs to 
an important family of tetrahedrally close-packed Frank- 
Kasper phases and is one of the closest low-order crys- 
talline approximants for dodecagonal quasicrystals which 
have similar local structural properties with the icosahe- 
dral glass (ICG) Q]- The latter implies that the study 
of the vibrational properties of the cr-phase should shed 
some light on similar properties in ICGs. Challenged 
by this possibility and motivated by a lack of available 
knowledge on dynamical properties of the real cr-phase we 
have carried out both, an experimental investigation as 
well as theoretical calculations, of the Fe-partial phonon 
density-of-states (DOS) for nominally Fe52.5Cr47.5 alloy 
having (a) a- and (b) cr-phase structure. 

The master alloy (ct-phase) was prepared by melting, 
in appropriate proportion, ^^Fe- enriched ('^95%) iron 
with chromium. The ingot was then cold-rolled down to 
a thickness of about 30/Km from which two 5x5 mm^ 
plates were cut out. They were next solution- treated 
at 1273 K for 72 h. One of the samples was afterwards 
transformed into the cr-phase by the isothermal anneal- 
ing at 973 K for 7 days. The verification of the transfor- 
mation into the cr-phase caused by such thermal proce- 
dure was done by recording a ^^Fe Mossbauer spectrum 
at 295 K. Experiments were conducted at sector 3-ID of 
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TABLE I: Atomic crystallographic positions and numbers of 
NN atoms for the five lattice sites of the Fe-Cr a-phase. 



Site 


Crvstalloeraphic positions 


NN 






A 


B 


C 


D 


E 


Total 


A 


2i (0, 0, ) 




4 




4 


4 


12 


B 


4f (0.4, 0.4, ) 


2 


1 


2 


4 


6 


15 


C 


8i (0.74, 0.66, ) 




1 


5 


4 


4 


14 


D 


8i (0.464, 0.131, ) 


1 


2 


4 


1 


4 


12 


E 


8j (0.183, 0.183, 0.252) 


1 


3 


4 


4 


2 


14 



the Advanced Photon Source. The vibrational properties 
such as the Fe-partial phonon density-of-states were stud- 
ied using the method of nuclear resonant inelastic X-ray 
scattering (NRIXS) Synchrotron radiation X-rays 

were monochromatized to a bandwidth of 1.2 meV and 
tuned in energy ranges of -I- /- 80 meV (room temperature 
measurements) and -20 meV to -1-80 meV (measurement 
at 20 K) around the ^^Fe nuclear transition energy of 
14.4125 keV. The X-ray flux and beam size at the sam- 
ple position were 4 x 10^ photons/s and 0.3x2 mm^, re- 
spectively. Data collection times were about 3 hours for 
the room temperature measurement of each sample and 
about 12 hours for the low temperature measurement of 
the cr-phase sample. We followed previously described 
evaluation proceduresjl, 0] using the publicly available 
PHOENIX software [7|. The following quantities were 
derived directly from the data: Lamb-Mossbauer factor, 
/, kinetic energy per atom, Ek, and the mean force con- 
stant, D. No specific assumptions about the character 
of the vibrations had to be made to obtain these values. 
The Fe-partial DOS was derived by direct data inversion 
using the Fourier-Log method under the assumption of 
quasi- harmonic vibrations. The consistency of this pro- 
cedure was verified by independent calculation of Lamb- 
Mossbauer factor, kinetic energy per atom, and mean 
force constant from the DOS and by agreement of these 
values with same quantities obtained directly from the 
data. Then the following quantities were calculated using 
the Fe-partial DOS: vibrational specific heat at constant 
volume, Cv, and vibrational entropy, S. The assignment 
of Debye temperatures, 0_d, is based on the Debye model, 
i.e., the DOS is proportional to energy squared, and they 
are widely used in the literature. With the determina- 
tion of the Fe-partial DOS, we have surpassed the Debye 
model but find it useful to provide Debye temperatures 
for comparison. Using the values of / and Cy, we deter- 
mined commonly presented values for Qd- 

In calculations, both phases of Fe-Cr alloy were mod- 
eled by the appropriate atomic configurations placed in 
a supercell with the periodic boundary conditions. The 
disordered a-Fe52.5Cr47.5 alloy was approximated by the 
a-FesoCrso one, for which we used the 2x2x2 bcc cell 
with 16 atoms. For random distribution of atoms, there 
are about 500 different atomic configurations to be con- 



sidered. However, for the sake of computer and time ca- 
pacity, we have chosen at random only five to be included 
in our calculations. The real cr-Fe52.5Cr47.5 sample was 
approximated by a cr-Fe53.3Cr46.7 one. The latter was 
studied in the 1x1x1 tetragonal supercell with 30 atoms 
(16 Fe and 14 Cr atoms). The structure optimization 
was done using the spin-polarized density functional to- 
tal energy calculations performed within the generalized 
gradient approximation (GGA) using the VASP pack- 
age 0, 0| • The valence electrons for each atom (electron 
configuration: dPs^ and d^s^ for the Cr and Fe atoms, 
respectively) are represented by plane wave expansions. 
The wave functions in the core region are evaluated us- 
ing the full-potential projector augmented-wave (PAW) 
method l3, 11 \ ■ The integrations in the reciprocal space 
were performed on the 8x8x8 and 4x4x4 grid for 
a- and tr-phase, respectively. During the optimization, 
the Hellmann-Feynman (H-F) forces and the stress ten- 
sor were calculated and the structure optimization was 
performed in two steps. First the lattice constants were 
determined assuming the appropriate symmetry, then the 
atomic positions were found in a fixed unit cell. The 
crystal structure optimization was finished when resid- 
ual forces were less than 10^^ eV/Aand stresses were 
less than 0.1 kbar. The calculated lattice constants are 
a = 5.64Aand a = 8.6OA, c = 4.74Afor the cubic and 
tetragonal symmetry, respectively. The optimized mag- 
netic moments on the Fe atoms are ordered ferromag- 
netically with the average values 2.03 and 0.96 /is, in 
the a and a phase, respectively. On the Cr atoms, the 
antiparallel arrangement occurs with the negative mean 
magnetic moments -0.17 and -0.28 in the respective 
phases. 

For the optimized structures the phonon dispersions 
and density of states were calculated using the direct 
method E^. The dynamical matrix of the crystal 
is constructed from the H-F forces generated while dis- 
placing atoms from their equilibrium positions. For con- 
sidered structures each atom must be displaced in three 
directions. For the a and cr-phases, a complete set of H- 
F forces is obtained from 48 and 90 independent atomic 
displacements, respectively. The amplitude of the dis- 
placements equals 0.03 A. To minimize systematic errors 
we applied displacements in positive and negative direc- 
tions. Finally, the phonon frequencies are obtained by the 
diagonalization of the dynamical matrix for each wave 
vector. The phonon DOS is calculated by the random 
sampling on the k-point grid in the reciprocal space, and 
then the thermodynamic functions are obtained within 
the harmonic approximation. 

The phonon DOS's measured on ^^Fe for both phases 
of the Fe-Cr alloy are presented in Figure 1. The dif- 
ferences in the energy range covered by the spectrum 
and the discrete structure of the spectrum are signifi- 
cant. The spectrum obtained for the a- phase is found to 
be similar to that of pure Fe [14| exibiting a distinct peak 
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FIG. 1: (Color online) Phonon DOS as measured on ^^Fe 
atoms at 298 K for the a- (circles) and the cr-phase (triangles) 
on the Fe52.5Cr47.5 samples. 



FIG. 2: (Color online) Phonon DOS as measured on ^^Fe 
atoms at 298 K (circles) and as calculated (dotted line) for 
the a-Fe52.5Cr47.5. 



at 36 meV. The Fe-partial DOS spectrum of the cr-phase 
demonstrates the additional high-frequency peak at 40 
meV, that is not observed neither in a-FeCr nor in pure 
bcc-Fe. There is also the shift downward of the low- 
energy peak, which makes the entire spectrum broader 
than in the a-phase. In Figure 2 the measured and calcu- 
lated phonon DOS spectra of a-FeCr are compared. The 
shape of both spectra is similar and two characteristic 
peaks of measured spectrum are reproduced satisfacto- 
rily. As in the chosen supercell, there are eight Fe atoms, 
hence by taking into account five different atomic config- 
urations, one considers vibrations of forty independent 
Fe atoms. Their partial DOS turned out to be different, 
but their shape was not correlated with the particular 
nearest-neighbor (NN) - next NN shell. Consequently, 
the final DOS was calculated using the partial contribu- 
tions with the same weights. The discrepencies between 
measured and calculated spectrum are likely caused by 
an incomplete representation of possible atomic configu- 
ration of FeCr disordered alloy in our model. 

Likewise, the calculations of cr-FeCr performed for only 
one configuration yield the phonon DOS exhibiting char- 
acteristic features of the experimental spectrum (Figure 
3). Observed discrepancies, like underestimation of the 
intensity of low-frequency peak or shift of high-frequency 
peak, are not significant. Using the theoretical result, it 
is feasible to separate the individual contributions to the 
total DOS generated by the Fe atom placed at each par- 
ticular crystallographic position. 

Thus, we see that the Fe atoms on sites A and C are 
causing the high-energy contributions to the DOS. The 
NN sites of these Fe atoms are placed at distances shorter 
than 2.48A. In pure bcc-Fe, all 8 NN atoms are situated 
at the same distance of 2.485A. The shorter distances 
between atoms result in the larger interatomic interation. 
Therefore, in cr-FeCr, the phonon frequencies higher than 
40 meV, which is a high-frequency limit of bcc-Fe DOS 



spectrum, are observed. Also due to larger dispersion of 
distances in the cr-phase, its phonon spectrum is much 
broader than in the a-phase, where the atomic positions 
are very close to those of the pure bcc-Fe. 




E [meV] 

FIG. 3: (Color online) Phonon DOS as measured on ^^Fe 
atoms at 298 K on the a- Fe52.5Cr47.5 sample (triangles) and 
as calculated (dotted line). DOS-curves for particular crys- 
tallographic sites are indicated, too. 

The data derived from the experiment and the theoret- 
ical calculations are displayed in Table I. One can see a 
very good agreement for some quantities like specific heat 
and quite good for other ones like force constant and en- 
tropy. The values of the Debye temperature derived from 
the specific heat are close to those calculated from the 
Lamb-Mossbauer factor, and there is a small difference 
for this quantity found for different phases. These find- 
ings disagree with experimentally found O^-values for 
the two phases using second-order Doppler shifts from 
Mossbauer spectroscopy [3]. In those studies, the Qd- 
values for the a-phase were larger than the ones for the 
cr-phase. This apparent discrepancy is explained by the 
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TABLE II: Physical quantities derived from the measured and 
calculated Fe-partial DOS. The units are meV/atom for vi- 
brational kinetic energy, Et; fcs/atom for vibrational entropy, 
S, and specific heat Cv', N/m for mean force constant, D; K 
for Debye temperatures, Od- 





experiment 


theory 


quantity 


a 


a 


a 


(J 


/@298K 


0.782±0.001 


0.768±0.002 






£;fc@298K 


42.4±0.1 


41.4±0.2 






D@298K 


156±1 


157±2 


149.41 


150.59 


Cv@298K 


2.747±0.006 


2.748±0.007 


2.746 


2.752 


S'@298K 


3.252±0.006 


3.347±0.007 


3.354 


3.412 


/@20K 




0.9150±0.0002 






£;fc@20K 




18.93±0.06 






L>@20K 




155.1±0.7 






/QOK 


0.9194±0.0001 


0.9149±0.0003 






EkmK 


19.32±0.07 


18.95±0.06 






eD(/@298K) 


417 


398 






eoiCv) 


399 


398 






ei3(/@20K) 




387 






ei5(/@0K) 


403 


385 







fact that the center shift measured by Mossbauer spec- 
troscopy constitutes a sum of the chemical isomer shift, 
which is independent on the atomic motion, and the 
second-order Doppler shift, which is a relativistic cor- 
rection to the atomic energy levels purely due to mo- 
tion. The second-order Doppler shift is proportional to 
the vibrational kinetic energy of the ^'^Fe atom [15|, i.e. 
SOD[mm/s] = -0.00565xi?fc[meV]. Only under the as- 
sumption that the chemical isomer shift is temperature 
independent, the Mossbauer measurement can provide 
the correct Debye temperature. This does not seem to be 
the case here, and our NRIXS data present a notable im- 
provement in the understanding of the role of vibrations 
in the Fe-Cr system. The value of the vibrational entropy 
measured at 298 K for the a-phase agrees quite well with 
inelastic neutron scattering results on samples of similar 
composition Also the difference in the entropy val- 
ues, AS = S„ - Sa = 0.095 ± O.OOgfcs as determined in 
the present experiment from the Fe-partial DOS agrees 
well with the corresponding difference calculated from 
the equation AS* = Sfcs ln(6Dcr/9DQ) = O.OTfcs, where 
ks is the Boltzmann constant and Qdi is the Debye tem- 
perature as determined from the second-order Doppler 
shift for the a {i = a) and the a {i = a) phase, respec- 
tively 0. The corresponding theoretical value of AS" is 
equal to 0.058 fc^. Taking into account the approxima- 
tions used in the calculations, the agreement seems to be 
quite satisfactory. 

In summary, we have revealed, both experimentally 
and theoretically, significant differences in the partial-Fe 
phonon DOS of the a and a phases of a quasi-equiatomic 



Fe-Cr alloy. Although, the sigma phase is a very complex 
object due to a high number of atoms per unit cell, five 
different sublattices with high coordination numbers (12- 
15), each showing chemical disorder, which altogether 
results in a huge number of possible atomic configura- 
tions, it was described reasonable well in terms of only 
one adequately chosen configuration. From the calcula- 
tions, it is also evident that the dynamics in particular 
sublattices is different. The method was also successfully 
used to calculate the dynamics of the disordered alloy in 
the alpha-phase. Although its crystallographic structure 
is much simpler, but, due to the chemical disorder, the 
number of possible atomic configurations is high, making 
the calculations not trivial. We have also obtained rele- 
vant thermodynamic quantities without necessity of us- 
ing empirical parameters. Such a complex alloy has been 
studied for the first time within the combined NRIXS 
and theoretical ab initio approach, and it may provide 
understanding of lattice dynamics in a wide variety of 
disordered systems. 
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